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Method for arbitrary phase transformation by a slab based on transformation optics 

and the principle of equal optical path 
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The optical path lengths travelled by rays across a wavefront essentially determine the resulting 
phase front irrespective of the shape of a medium according to the principle of equal optical path. 
Thereupon we propose a method for the transformation between two arbitrary wavefronts by a slab, 
i.e. the profile of the spatial separation between the two wavefronts is taken to be transformed to a 
plane surface. Interestingly, for the mutual conversion between planar and curved wavefronts, the 
method reduce to an inverse transformation method in which it is the reversed shape of the desired 
wavefront that is converted to a planar one. As an application, three kinds of phase transformation 
are realized and it is found that the transformation on phase is able to realize some important 
properties such as phase reversal or compensation, focusing, and expanding or compressing beams, 
which are further confirmed by numerical simulations. The slab can be applied to realizing compact 
electromagnetic devices for which the values of the refractive index or the permittivity and perme- 
ability can be high or low, positive or negative, or near zero, depending on the choice of coordinate 
transformations. 

PACS numbers: 42. 15. Dp, 41.20.Jb, 42.79.-e, 02.40.-k 
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INTRODUCTION 



In the last few years, transformation optics [l|,l2[, based 
on the form- invar iance of Maxwell's equations under co- 
ordinate transformations, has been attracting more and 
more attention. According to the equivalence of space 
geometry and material in routing light, it provides a ro- 
bust method to control electromagnetic waves by ma- 
terials usually implemented by metamaterials This 
has led to remarkably successful applications to design- 
ing electromagnetic devices, which may be mainly clas- 
sified into four classes. The first is to control the path 
of light ray, such as invisibility cloaks 0-|6j|, light con- 
centrator/absorber @J8[, beam expanders/ compressors 
HJ GU, beam shifter [lj [H, and beam bends [H, Q. 
The second is to manipulate the polarization [l5| and the 
third is to control the amplitude [l6l-[l8j. While the last 
is to mold the phase with notable exam ples as the field 
rotators [19| and phase transformers [2Ch22(. 

There exists a limitation on the wave through most of 
the above devices yet. That is, the exit wave would return 
0] or be parallel to the incident wave direction, or be 
perpendicular to the exit interface [13j because the wave 
vector or phase front has not been changed by transfor- 
mations. Recently, phase transformers in the last class 
achieved the direction of exit wave different from the in- 
coming one, but at the cost of irregular profiles due to 
the direct transformation used. It is natural to ask if a 
device can be designed independent of the geometrical 
shape, such as a slab, to manipulate the exit wavefront. 
As well known the planar structure is preferred in practi- 



cal applications, especially as plug-and-play devices, be- 
ing compact enough as well as free of aberrations (23| . 
Indeed certain planar designs have been used to realize 
directive emission with high performanc e by transform- 
ing cylindrical waves into plane ones |24j29j. However a 
general theory for phase transformation through planar 
configurations has not been proposed yet. 

In this work we introduce a new coordinate transfor- 
mation method from the fundamental concepts of wave- 
front and optical path length (OPL) in ray optics to de- 
sign a slab to manipulate the phase. Through the slab 
the conversion between two arbitrary wavefronts can be 
achieved. The slab can be realized by materials for which 
the refractive index or the permittivity and permeability 
can be high or low, positive or negative, or near zero, de- 
pending on the choice of transformations. The results can 
be applied to realizing compact electromagnetic devices, 
such as wave deflector, flat lens, phase compensator and 
beam expander or compressor. In addition, the method 
takes a straightforward form and can completely avoid 
possible singular points in constitutional parameters. 

This work is organized as follows. In Sec. Ill A I we de- 
rive from the concepts of wavefront and OPL the coordi- 
nate transformation method related to planar wavefronts, 
i.e. the inverse transformation method, then extend it to 
a general form in Sec. Ill Bl In Sec. IIIII we realize three 
kinds of phase transformation as an application of the ap- 
proach and discuss the influence of the choice of coordi- 
nate transformation on material parameters. Discussion 
and conclusion are made in the last section. 



II. THE BASIC FORMULATION 
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Let us begin with the concept of wavefront and OPL. 
As well known the wavefront is the surface of constant 
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phase and the phase differences of any two correspond- 
ing points on two sequential wavefronts are identical [30| . 
The optical path lengths (OPLs) involved are also equal, 
stated by the principle of equal optical path [3l|. In ac- 
cordance with Fermat's principle, the rays across initial 
wavefront traverse individual paths so that the OPLs or 
the times of travel are extrema, whether through conven- 
tional optical components [32[ or newly developed trans- 
formation media. If these OPLs are distinct, they arrive 
the exit surface at different moments and the associated 
phases are varied. Thereby the resultant wavefront is 
distorted according to the principle of equal optical path. 
So, strictly speaking, it is the OPLs travelled rather than 
the shape of medium that fundamentally determine the 
resulting phase front (3l|, [32|. Theoretically, therefore, 
one phase transformation can be realized by materials of 
various profiles as long as the needed OPLs, i.e. phase 
differences, are satisfied. On the other hand, the trans- 
formation between wavefronts corresponds to one specific 
coordinate transformation in a mathematical perspective 
[l9l-[29|. Thus it is entirely feasible to construct ma- 
terial with a slab configuration to transform wavefront 
by transformation optics. These form the basis for our 
method. 



A. The inverse transformation method for the 
transformation associated with planar wavefronts 

For simplicity we consider the two-dimensional wave 
problem and first establish the framework of transfor- 
mation to convert a planar wave front into any desired 
phase front. To achieve that the usual approach in lit- 
erature is to transform a plane surface to the desired 
wave shape immediately, i.e. from B'C to A'E' as 
shown in Fig. DJa), which may be called direct trans- 
formation method (DTM) [1h22|. Therein the virtual 
region (OB' CD, enclosed by dotted lines) is transformed 
into the physical one (OA'E'D in yellow). Such a DTM 
is intuitive enough but results in designs with irregular 
shapes, OA' E'D. Instead, we adopt an inverse transfor- 
mation method (ITM): the reversed shape of the desired 
wavefront is converted to a planar one, i.e., from AE to 
BC, whereby the design is a slab, OBCD, as shown in 
Fig. DJb). Here it is the virtual region OAED that is 
transformed into the flat physical one OBCD. 

In the transformation of any point from B'C to A'E' , 
the distance travelled is g(y') + a' — b' , where B'C is set 
as x = b' and A'E' is defined as x' = a' + g(y'). Instead, 
in the transformation way from AE corresponding to x = 
a — g(y) into BC with x = 6, the separation is g(y) + b — 
a. Here (x, y) is an arbitrary point in the virtual space 
while (x',y') in the physical space. Practically constants 
a' — b' or b — a do not affect the resulting shape of phase 
front, so the effective OPLs are identical between the two 
ways since y = y' . Thereupon the ITM can implement 
exactly the same function as DTM. In detail, the slab in 
Fig. QJb) can convert a planar wave front into a curved 
one, the opposite shape of AE, which is just the function 
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FIG. 1: Schematic diagram of the space transformation to 
convert a planar phase front into a curved one by (a) direct 
transformation method in literature and (b) inverse trans- 
formation method adopted in the present work. AE is a 
reversed version of A'E'. The dotted regions are the vir- 
tual spaces while the real-lined regions indicate the physical 
spaces. OA' = a', OB' = b',OA = a, OB = b, and OD = d. 



of OA'E'D in Fig. QJa). Meanwhile, within the context 
of DTM the slab should be capable of converting a curved 
wave AE into a plane wave BC. Therefore the slab by 
ITM has a twofold function which we denote by I and 
77, respectively. For simplicity, we do not consider the 
transformation in the y direction. 

By ITM the spatial distortion can be recorded as a 
mapping between the original and transformed spaces, 



x' p(x,y) 



bx 
A' 



z = z, 



where p(x, y) is a scaling factor and 
A = a-g(y) 



(1) 



(2) 



corresponding to AE, the reversed shape of the wanted 
wavefront. The mapping means that the virtual space 
is stretched or squeezed proportionally into the physical 
space. Theoretically one can choose a transform freely 
[HQ. The present choice is made as p(x,y) = 1 in order 
that the parameters change linearly and the slab can be 
implemented easily. The difference from the literature 
is that it has such a simplicity that it needs not con- 
sider the details of transformation. According to trans- 
formation optics, the permittivity tensor e and perme- 
ability tensor \i in the transformed coordinate system 
are connected with the original e and ti Q by the rela- 
tionships: e — Je Q J T /det(J), fj, = J/j, J T /det(J) where 
J\ = dx % I dx % is the Jacobian transformation matrix [3[ . 
Applying Eq. (pQ) gives a general result of the relative ma- 
terial parameters for the conversion between plane wave 
and arbitrarily curved wavefront: 
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where A' = dA/dy and the original space is considered 
as vacuum. Notice that the superscripts of coordinate 
variables have been omitted for simplicity. No singularity 
exists for e and as long as AE does not intersect with 
?/-axis. Thereupon, the problem of singularity [2?], HH 
arising mathematically can be avoided completely by tak- 
ing appropriate transformation ways. It is important to 
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FIG. 2: Schematic diagram of the space transformation to 
convert a curved phase front into another curved one by (a) 
direct transformation method in literature and (b) the trans- 
formation method in the present work. AE, a profile of the 
spatial separation of the two wavefronts B'C' and A'E' , is 
taken to be transformed into a plane. 



note that the complex amplitude transmittance [33| of 
electromagnetic wave through the slab can be written as 



t 



-iko A 



(4) 



where fco is the free-space wavenumber and the constant 
phase factor has been omitted. 



B. The general transformation method for the 
transformation between two arbitrary wavefronts 

In an analogous fashion to the above discussion, we 
now examine how to convert a curved wave front into 
any desired one by a slab. Fig. [2] applies as a sketch 
wherein B'C is curved now. The DTM in literature is 
to transform a curved wave to the desired wave shape im- 
mediately, i.e. from B'C to A! E' , and results in designs 
with irregular shapes, OA' E'D. Here, we generalize the 
ITM following the same principle based on OPL. Briefly, 
one also can imagine compressing A'E' into a plane BC, 
then the original B'C would undergo the same distortion 
in the x direction becoming AE. However, the transfor- 
mation becomes that the profile of the spatial separation 
between the original and desired wavefronts is converted 
to a plane surface, i.e. from AE to BC, to realize the 
design by a slab as OBCD in Fig. Efb). 

In transforming any point from B'C to A'E', the 
distance experienced is #2(2/) — 9i{y') + a' — b' , where 
B'C and A'E' are supposed to be x' = b' + gi(y') and 
x = a' + #2(2/), respectively. Instead, in the transforma- 
tion way from AE with x = a + g\{y) —#2(2/) to BC with 
x = 6, the separation traversed is gi{y) — #2(2/) + b — a. 
Actually a' — b' or b — a do not affect the resulting shape 
of front, so the effective OPLs are identical for the two 
ways since y = y' . Therefore the results Eqs. (pQ) and ([3]) 
still hold true except that A now becomes 



a + 9i(y) -92(y) 



(5) 



corresponding to the profile of spatial separation between 
the two wavefronts. If the incident wave front is planar, 
e.g. a + gi(y) = b', then A = b' — #2(2/), which is the 
reversed shape of the outgoing wavefront. That is to say, 
to transform a plane wave into a curved wavefront by a 



slab, it only needs to transform the reversed shape of the 
latter into the former. This is just the above conclusion 
of Eq. (j2j), so Eq. (j5]) is a general result. 

It is important to note that in deriving Eq. (|2]) or 
Eq. (j5j), AE should be shifted so as not to intersect the y 
axis to avoid the singularity in e and fj, in Eq. (|3]). That 
is, A can be added with a constant that will not change 
the generated wavefront, but only affect the values of e 
and 11 in the slab. We will address this point later. 



III. APPLICATION OF THE 
TRANSFORMATION METHOD 

Applying the above method, a wavefront can be con- 
verted into any desired one through a slab with e and \i 
determined by Eq. (|3j). In the following we discuss three 
kinds of phase transformation to illustrate the method. 
Incidentally, similar wavefronts could have been gener- 
ated by other media using different methods somewhere 
in literature, e.g. conventional lens or transformation 
media, but our principal emphasis is on the new general 
method used to design flat media to induce phase shift 
and bend wavefront. 



A. Planar to planar wavefront transformation. 

Suppose a plane wave to be deflected into other direc- 
tions by the slab. Without loss of generality, the reversed 
wavefront AE is chosen as y — — (x — a) cot 9 such that 
the exit wave is deflected clockwise by an amount 0. Then 



A = a — y tan 6 



(6) 



and using Eq. (J3j) one obtains the relative material pa- 
rameters for the slab. 

To test the effectiveness of the ITM and the proposed 
design, suppose a = 0.4m, b = 0.3m and = ir/6. Con- 
sider a z-polarized Gaussian beam with a wavelength of 
A = 0.05m incident on the slab of d = 0.8m (the choices 
are uniform throughout), the electric field distribution 
obtained from a finite-element simulation is shown in 
Figure [3] We see that the incident wavefront is grad- 
ually tilted inside the slab until the deflection result is 
transferred out to the free space. The exit beam acquires 
a new direction of propagation unlike that of the inci- 
dent one and the deflected wavefront is not necessarily 
parallel to the exit interface. Interestingly, by comparing 
(a) and (b) we find that the beam is deflected by the 
same amount whether for normal or oblique incidence. 
The underlying mechanism is that this transformation is 
in fact an operation of rotation on the local coordinate 
system which will be imposed on any fixed incident field 
vectors, resulting in wave vectors with identical deflec- 
tion. This fact confirms the twofold function of the slab: 
/-Convert a normally incident plane wave into another 
with the deflected angle //-Convert an oblique plane 
wave with the incident angle — into a normal one. Also 
the ITM is justified. 
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FIG. 3: Normalized i^-field patterns in the slab deflector for 
(a) normal incidence with exit angle ty/6 and (b) incident 
angle — n/6 with normal outgoing, and (c) in a splitter formed 
by the upper slab with deflection angle of — 7r/12 and the lower 
of 7r/6. (d) The permittivity at the center of the slab varies 
with the width b. 




FIG. 4: Normalized i^-field patterns in the slab applying (a) 
converging / '(b,c,d) diverging surface to planar one transfor- 
mation. The slab can convert a planar wavefront to a (a) 
diverging /(b,c,d) converging one. The surfaces to transform 
are chosen as A = 0.05 + (y - 0.5) 2 in (a), 0.2 - (y - 0.5) 2 
in (b), -0.2 - (y - 0.5) 2 in (c), and -(y - 0.5) 2 in (d). The 
incident beam centers are (0, 0.5) and the waists wq are 0.2m. 



To change the deflection amount, one can adjust the 
parameter 0. Positive or negative i.e. clockwise 
or counter-clockwise deflection, even extreme angles ap- 
proaching the boundary, can be realized. Alternatively, 
one can realize different deflections by stacking certain 
pieces of slabs with given deflection amounts. In addi- 
tion, combine two such slabs and a beam splitter is cre- 
ated, e.g. Fig. [3jc). A benefit of this splitter is that 
a beam can be split into any two directions, provided 
that the deflection angles are set accordingly. We also 
find that the width b can be chosen freely, so the slab 
can be made as compact as desired only at the expense 
of a larger range of material parameters, as shown in 
Fig. Bid). 



B. Planar to curved wavefront transformation. 

For example, we consider how to convert a plane wave 
to a diverging/converging wave. AE is assumed to be a 
paraboloid 

A = a-y 2 /2R(a) (7) 

positioned at x = a where the radius of curvature 

R(x) = ±x(l + x 2 /x 2 ) (8) 

and the Rayleigh range xq — w^tt/X corresponding to 
a Gaussian beam of waist radius wq. Eqs. Q and (|8j) 
will be of continuous use from now on. Applying the 
ITM, the inverse transformation is taken from a converg- 
ing/diverging surface to a planar one. Then the material 
parameters are obtained by using Eq. (|3j). Similar trans- 
formations were used to generate plane waves [25|, l2rjj , 
but here we make use of the map inversely. 

In order to show the phase evolution clearly, we choose 
a — 0.05m/0.2m. The other values are wq = 0.2m, 



R = — 0.5m/0.5m and the waist center (0, 0.5) for the in- 
cident beam which approximates a planar wavefront, and 
the simulated results are shown in Fig. |4] (a) /(b). Note 
that we have adopted the convention that a converging 
wavefront has a negative radius of curvature whereas a 
diverging wavefront has a positive radius of curvature. In 
the slab one can see the incident planar wavefront grad- 
ually evolves into a diverging/ converging one at x = 
i.e., the function /. The central area of the wavefront 
travels faster /more slowly than its outer edges and then 
the midregion overtakes /lags behind those edges, thereby 
continuously bending the wavefront. In terms of Fer- 
mat's principle, it means that the refractive index is 
lower/larger in the central than in outer edges. In Fig. [T] 
the generated surface A!E' is the mirror image of the 
transformed surface AE, so one can see that the radius 
of exit phase surface is 0.5m/— 0.5m just contrary to R 
and the new waist varies accordingly. 

On the other hand, the slab in (a) / (b) enables II- 
make a converging/diverging wavefront bend into a pla- 
nar one, as shown in the literature applying similar coor- 
dinate transformations [2^-[27(] . Hence the slab also has a 
twofold function which may find import applications. For 
instance, it is capable of implementing fiat lens. Recall 
that for a thin lens the complex amplitude transmittance 
is t = e tk onbo-ik y /2f ^ w h ere n i s the refractive index, bo 
is the thickness of the center, and / is the focal length 
[33( | . Therefore, the slab determined by Eq. ([7]) is equiv- 
alent to a thin lens with the focal length f = R, which 
was confirmed by Fig. 2J Apart from its shape being a 
compact slab, such a lens has another advantage that the 
location of focus can be tuned continuously by changing 
the material parameters or discretely by stacking cer- 
tain pieces of slab with unit focus distances. Note that 
the result in (b) is very similar to but nevertheless dis- 
tinct from previous results of flattening lens [34|, [35| . The 
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transformation in the latter was to reshape a known opti- 
cal component, whereas the present method immediately 
take the wavefront to transform and thus is endowed with 
great flexibility, so results e and \i are different, too. 

To examine the role the factor a plays in the result, 
we alter its value in (b) as a — —0.2m and 0m, and the 
results are shown in (c) and (d), respectively. Note that 
the original wavefront in (c) and (d) are chosen left to the 
slab, whereas that in (b) is inside the slab. Comparing 
(b) and (c) or (d), one can readily see that the period 
number of wave inside the slab equals a/ A. This fact can 
be understood by Fig. [1] where the original wave distri- 
bution in OAED is transformed into that in OBCD and 
OA = a, so the period number of the original wave along 
OA is just that in the slab along OB. Particularly in- 
teresting is the crest in (d) that becomes quite elongated 
and stretches into a long beam inside the slab, exhibit- 
ing a extremely large wavelength. Correspondingly, the 
index of refraction should be nearly zero [36| , which will 
be proved later. 

At the same time, it can be seen from (b), (c) and (d) 
that the exit wave fronts remain fixed, converging at a 
point R away from the slab. This fact results from that 
the curvature of radius in the transform is uniform. So 
they can implement the same function and thus the orig- 
inal surface AE in Fig. [I] need not be chosen in the slab. 
Particularly, in contrast with (b), it seems like that the 
central region of wavefront travelled faster than the outer 
edges inside the slab in (c) and (d), like in (a). However, 
the fact is quite the contrary and the true reason is that 
negative refraction occurs in (c) and (d), thereby intro- 
ducing negative phase shift. Since the central area of 
the beam travels more slowly than its outer edges, the 
subtracted phase is less and the accumulated phase of 
the midregion is advanced with respect to those edges 
on the back surface, whereby the outgoing wavefront are 
deformed as in (b). This point is evident from the oppo- 
site phases traversing the boundary in (c) and (d). Such 
a process leads to the reversal of phase that may find 
important applications, such as wavefront correction in 
imaging or adaptive optics. In contrast the phases across 
the interfaces are uniform in (a) / (b) and the refractive 
indexes are positive. Therefore, we come to an impor- 
tant conclusion that the ways to get a desired wave are 
diverse and then the same to the choices of slab, thereby 
increasing the possibility of realization. 



C. Curved to curved wavefront transformation. 

To be specific we investigate the conversion of a diverg- 
ing wave x = a + gi(y) with a = 0.3m, R = 0.5m and 
the beam center (—0.15,0.5) to a converging one #2(2/) 
with R' = —1.2m, where Eqs. ([7j) and © have been ap- 
plied. Then the material parameters are obtained by us- 
ing Eqs. (j3j) and (j5j), and the simulated results are shown 
in Figs.[5fa) and (b). Apparently the left diverging wave 
traversing the slab is deformed into a converging one. 
As in the above subsection, the period number of wave 




FIG. 5: Normalized i^-field patterns (a) and the energy den- 
sity (b,c,d) in the slab transforming a diverging wavefront to 
a converging one. (c) displays the focus on the exit boundary 
when —a equals to the distance of beam center from the slab. 
If —a becomes larger, then another focus will appear inside the 
slab, leading to focusing of light in (d). The two wavefronts 
are chosen as (a,b) a + gi(y) = 0.3— (2/ — 0.5) 2 , (c) a + gi(y) = 
-0.15 + (y - 0.5) 2 , (d) a + gi (y) = -0.3 + (y - 0.5) 2 , while 
#2(2/) = (y — 0.5) 2 /2.4 and incident beam center (—0.15,0.5). 



inside the slab equals a/X and R' determines the gener- 
ated wavefront radius of curvature. At the same time, the 
slab induces a transmittance the same as that of a thin 
lens with the focal length satisfying 1/ f = 1/R — 1/R' '. 
Moreover, the focus can be tuned continuously by chang- 
ing the parameters of transformation. For example, let 
—a equals to the distance of beam center from the slab, 
then the beam will focus at the exit surface, as shown 
in Fig. [5fc). The reason is that the choice of transfor- 
mation is to let the wavefront at —a, i.e. incident beam 
center, be converted to a plane surface. The slab now 
acts as a phase compensator that reverses the amount 
of phase change within the region a away from the slab. 
Therefore, it is expectable that, if —a is larger than the 
distance of beam center from the slab, then the beam 
will focus inside the slab, as shown in Fig. [5fd). Such 
a function can be utilized as beam relaying and flat lens 

As known beam expansion or compression has been 
achieved by manipulating the longitudinal energy flow 
[HUH; that is, expanding or squeezing the virtual space 
in the y direction in Fig.[TJ In contrast, the same task can 
be performed now by the transformation on the trans- 
verse phase, that is, distorting the virtual space in the 
x direction. From the above, the wavefront with the ra- 
dius of curvature R at a is transformed into another one 
with the radius of curvature R' at b. The widths of the 
two waves, w and w' across the back interface can be 
considered unchanged. Then, for the outgoing wave, in 
accordance with R' and w' ', the beam center x away from 
the slab and the waist radius w' can be found, i.e., x = 
-R'/[l + (Afl'/W 2 ) 2 ] and w' = w' 2 /[l + (irvP/XR') 2 ]. 
Therefore, the width of the resultant beam can be engi- 
neered through choosing appropriate wavefront to trans- 
form and R' . In Figs. [SJa) and (b), the wavefront trans- 
formed with a distance from the center x — 0.45m has a 
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width w = wq<\/1 + (x/xo) 2 = 0.15m, so the waist radius 
of the resulted beam equals to w f = 0.1m, two times the 
original one. Evidently, the incident beam is expanded 
into another one with wider waist. The contrary condi- 
tion, beam compression, can be realized similarly. 



D. Influence of the coordinate transform on 
material parameters. 

Next we discuss the effect of the choice of transforma- 
tion on the material parameters. As the above analysis 
indicates, e.g. Fig. Efd), the values of e and fj, change 
as the slab thickness varies. On the other hand, the sign 
and magnitude of A also plays an important role in de- 
termining e and /J,. We have known from Sec. [IT] that to 
avoid the singularity in e or /L£, AE, i.e. A, can be shifted 
some amount along the x axis. Despite no impact on the 
generated wavefront, it undoubtedly influences the val- 
ues of e and /jl in the slab. From Eq. (|3j) we conclude 
that: (i) The larger the magnitude of A, the larger the 
range of values of e and fi. In terms of the principle of 
the equality of OPL, wave need to travel longer path A 
to generate the phase difference and then a width-fixed 
slab has to increase the magnitudes of e and /j,. (ii) If A 
becomes negative, e and \i will also take negative values. 
Especially, if A becomes zero, e and fi would approach 
infinity. In this case, A can be added a positive constant. 
As an example, we present the constitutional parameters 
in Figs. H3(a)-(c) corresponding to the slab in Figs.HJb). 

Further, in order to characterize the slab more clearly 
and comm only we introduce the effective average refrac- 
tive index (32, [33|, [3?], [HI in terms of the complex ampli- 
tude transmittance Eq. (J2J), 



eff 



A 
~b' 



(9) 



It indicates that the anisotropy is reduced such that non- 
diagonal elements in J are eliminated by choosing orthog- 
onal grids in virtual space to transform. Then the slab 
can be considered as nonmagnetic and be realized by di- 
electrics only though the performance may be sacrificed 
a bit [37], HI] . The effective refractive index on the trans- 
verse of the slab for different choices of a are shown in 
Fig.[6jd) which has included the results in Figs.HJb), (c) 
and (d). It shows that n e ff increases with a and that 
n e f f in the center region is larger than in the outer area. 
This is reasonable because the larger a, the more origi- 
nal space will be squeezed into the slab in Fig. QJb). For 
the midregion, if a < 0, then n e ff < 0, e.g., the left 
two lines in Fig. [6fd) corresponding to n e ff for slabs in 
Figs. 11(d) and (c), respectively. Accordingly, the choice 
of transformations determines whether the values of n e f f 
or e and /j, are high or low, positive or negative, or near 
zero. Therefore, it needs to take all factors into consider- 
ation and select appropriate values to realize the design 
by materials available. 
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FIG. 6: (a), (b) and (c) respectively corresponds to s xx (lixx), 
£yy(l^yy)/£zz(l^zz) and £xy(/J<xy)- (d) shows the effective re- 
fractive index on the transverse of the slab for different choices 
of a with b = 0.3. The three dashed lines in (d) from left to 
right correspond to n e ff for slabs in Figs. H^c), (d), and (b) 
respectively. 



IV. DISCUSSION AND CONCLUSION 

Before concluding this work it merits some further dis- 
cussions. First, to implement the above design, metama- 
terial is competent without any technological difficulties 
[39[ [iof . Now another simpler way is to further discretize 
the material parameters and build by dielectrics, which 
has been proved efficient to realize lossless and broad- 
band transformation devices 0, H, [HI, OS, El IH • Sec- 
ond, according to the principle of equal optical path the 
method can be generalized to realize arbitrary phase con- 
version through media with arbitrary geometric shapes. 
The coordinate transformation should guarantee that the 
needed OPLs are satisfied. Last, we only considered the 
transformation in two dimensions. Actually, the method 
is applicable to three-dimensional problems and more 
new phenomena may be found. For example, the method 
can be used to transform planar wavefront of Gaussian 
beam into helical one of vortex beam carrying orbital 
angular momentum [431 ]. 

In summary, we have proposed a planar phase trans- 
former to achieve the conversion between any two wave- 
fronts. From the point of view of OPL we present a gen- 
eral method of transformation wherein the profile deter- 
mined by the difference of the two wavefronts is taken to 
be converted to a plane surface. For the mutual conver- 
sion between planar and curved wavefronts, the method 
is greatly simplified as the ITM in which it is the oppo- 
site shape of the desired wavefront that is converted to a 
plane. Applying the method, three types of phase trans- 
formation are investigated, which can be used for wave 
deflector and flat lens, and are further confirmed by nu- 
merical simulations. In addition, we find that the trans- 
formation on phase can realize some important properties 
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such as phase reversal or compensation, focusing, and ex- 
panding or compressing beams. The slab can be realized 
by materials whose values of refractive index or e and \i 
depend on the choice of transformations. The method 
provides a new insight into transformation optics and we 
expect it will be of practical importance for designing 
compact optical devices. 
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